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Abstract 

O . 

, We have investigated the evolution of the neutrino asymmetry in active-sterile neu- 

^"*^ ' trino oscillations in the early universe. We find that there are large regions of parameter 

1 , space where the asymmetry is extremely sensitive to variations in the initial asymme- 

' try as well as the external parameters (the mass difference and the mixing angle). In 

these regions the system undergoes chaotic transitions; however, the system is never 
C*") ' truly chaotic in the sense that all information about initial conditions is lost. In some 

cases though, enough information is lost that the final sign of the neutrino asymmetry 
is stochastic. We discuss the implications of our findings for Big Bang nucleosynthesis 
(BBN) and the cosmic microwave background (CMB). 
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Neutrino oscillations have been proven beyond reasonable doubt by the observation of 
the solar and atmospheric neutrino anomalies. The solar neutrino problem can be explained 
by oscillations between v e and [J!], whereas the atmospheric neutrino deficit is quite 
\ nicely explained by oscillations between and v T ||. However, the LSND experiment also 

i ' claims detection of neutrino oscillations between and v e , but with a much larger mass 

difference than found from the solar neutrino experiments |J. There is no possible three- 
r-j ! neutrino solution to the combined data. Either one of the interpretations is wrong, or there 

is a fourth neutrino species, responsible for either the atmospheric or the solar neutrino 
anomaly. It has proven difficult to rule out the existence of such a fourth neutrino. It 
cannot interact via the usual weak interactions, or it would have been detected in Z-decay 
experiments. However, such additional light sterile neutrinos are predicted to exist in many 
extensions of the standard model, and it is certainly worthwhile to study the implications 
of such neutrinos for cosmology. 

Active-sterile neutrino oscillations in the early universe have been intensely studied for 
more than a decade. The pioneering studies concentrated on effects from excitation of inert 
degrees of freedom from the background plasma on BBN predictions of the abundances of 
the light elements [4-10]. Extra degrees of freedom exited before the BBN epoch increases 
the relativistic energy density of the universe. This modifies the expansion rate of the 
universe and changes the outcome of BBN, most notably the primordial helium abundance, 
Y P . 

It was realized early on that neutrino oscillations in a CP asymmetric plasma like the 
early universe could amplify an initial neutrino asymmetry (L u = (N u — Np)/N^) ||. 
Initially, the effect was thought unimportant m, and it was not until 1996 where Foot, 



Thomson and Volkas [11] showed that for for some mixing parameters, an important and non 
negligible asymmetry can actually be generated and affect BBN-predictions JTT] , |i~2fl . These 
results were confirmed by Shi [13], who in addition argued that in parts of the (5m 2 , sin 29)- 



parameter space, the sign of the asymmetry can go though a period of rapid oscillations 
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before settling down at a final value. The magnitude of the final asymmetry seems quite 
robust, but the final sign is very sensitive to small changes in the initial asymmetry and the 
mixing parameters. These results have later been investigated in further detail. 



Based on the quantum rate equations (QRE's), Enqvist et al. [14| have confirmed Shi's 
results The set of QRE's is an approximation to the full quantum kinetic equations (QKE's) 
for the system, based on the assumption that all modes have the average momentum, p = 
(p) ~ 3.15T. This approximation is excellent outside resonance regions, but its applicability 
inside resonances has been questioned. The reason is that in the QKE's different momentum 
modes pass through the resonance at different temperatures. Assuming only one mode could 
possibly lead to different behaviour of the oscillations inside the resonance. 

Unfortunately it is not yet possible to solve the full QKE's through a strong resonance 
due to numerical limitations [15, [l6|]. This fact has lead to recent controversy^] about the 



behaviour of the asymmetry when the full set of QKE's is employed. We have chosen 
to work with the quantum rate equations in the present treatment because the equations 
are numerically solvable and reliable conclusions can be drawn from them. However, we 
caution that some of our conclusions do not necessarily apply to the full set of quantum 
kinetic equations. 

The question of chaoticity is important since chaotic generation of lepton number might 
lead to the formation of leptonic domains. Not causally connected regions could develop 
different signs on the leptonic asymmetry. Diffusion of neutrinos over leptonic boundaries 
would be expected to suppress asymmetry and might also open a new channel of producing 
sterile neutrinos via resonant MSW conversion of active neutrinos at domain boundaries 

Our goal is to investigate to what extend the final sign of the asymmetry can be stochas- 
tic. In order to do this we employ a Lyapunov-exponent analysis of the system and calculate 
the loss of information about initial conditions experienced as the system passes through 
the resonance. Finally, we discuss the implications of our findings for BBN and the CMBR. 

Equations — The two- neutrino mixing is characterised by the vacuum mixing angle 9. 
The matter eigenstates are connected to the flavor eigenstates u a>s via the transforma- 
tions^] 

v a = cos 9 v\ — sin 9 v<i 

v s = sin 9 v\ + cos 9 v%. (1) 

If the two mass eigenstates have masses m\ and 777.2, the squared mass difference is defined 
as; bvn 2 = — mf. The objects of interest are the reduced 2x2-density matrices p v and 
Py for the neutrino and anti-neutrino ensembles respectively. Expanded in jj], ^(J Pauli spin 
matrices these take the form 

#, = ijb(l + P-(T); /7 77 =ipo(l + P^). (2) 

The coefficients in the expansion form the components of the polarisation vector P. 



3 The controversy raised by Dolgov et al. |17|] seems to have been resolved and the conclusions of Dolgov 
et al. rejected (see Refs. ]l8|). 

4 Note, the absolute phase of the non-diagonal terms is arbitrary; our choice differs from some in the 
literature. 
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The diagonal elements in p u are the relative number densities of the active neutrino and 
the sterile neutrinos normalised to their equilibrium values 

nu a = \Po(l + P t ), n Vs = ip (l - Pz) (3) 

with similar equations for the antineutrino ensemble. With the adopted normalisation, Pq 
then equals the sum of the relative number densities of the mixed neutrinos (2 in chemical 
equilibrium). P z measures the distribution between the active and the sterile neutrinos. If 
P z =l all neutrinos in the mixed ensemble will be of the active type, whereas all neutrinos 
will be of the sterile type if P z = — 1. 

Instead of solving the full set of QKE's for p u {p) and pp(p), we solve the quantum rate 
equations for p v (j> = (p)) and p v {jp = (p)), where (p) ~ 3.15T []. 

We focus on the part of parameter space where the the final sign of the leptonic asym- 
metry appears extremely sensitive towards small changes in the parameters. In this part 
of parameter space collision terms can be neglected so that the resonance is non-adiabatic 
and no significant amount of sterile neutrinos are produced. Po is then a constant which 
can be put equal to unity. The approximation is valid for mixing parameters obeying the 
inequality Q 

|£m 2 |sin 4 20 < 1(T 9 . (4) 
The coupled equations of motion then are 

P = VxP-D P T , F = VxP-D P T , (5) 

where P = c?P/dt, and the transversal polarisation vector Pt = P x x + P y y has been 
introduced. For definiteness we shall here focus on v T — v s oscillations; other cases are 
obtained from this by simple redefinitions. The damping coefficient D signifies quantum 
damping introduced in the system by coherent and incoherent scattering of neutrinos on 
particles in the background plasma. In the case of v T — u s oscillations it is given by ppfl 

D~D~ 1.83 G 2 F T 5 ~ 2.49 x 1(T 46 T 5 eV~ 4 . (6) 

The rotation or potential vector V arises from an expansion of the total weak Hamilto- 
nian in Pauli spin matrices. In the case of v T — u s oscillations it is given by 

V = -b s x+(-b c -V, T )z, V = -b s x+(-b c -V^)z, (7) 

where 



sin 29 5m 2 sin 29 5m 2 

s = 2(pj ~ 6.3 T ' 

and 

cos 29 5m 2 cos 29 5m 2 

Or = r- ; — . 

2{p) 6.3 T 



(8) 



(9) 



5 One can also argue that the momentum should be p ~ 2. IT, which is the maximum for the Fermi-Dirac 
distribution function. Taking this value instead does not make any qualitative difference, but can make some 
quantitative difference because the resonance temperature is altered by a factor (3. 15/2. f ) 1//2 ~ 1.22. 
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The term V„ T (V Ut ) is an effective energy contribution to the i> T {y T ^) from neutral current 



and charged current interactions with particles in the background. It is given by [21 



V T = -br + b A 5P z , V T = -b T - b A 5P z (10) 

The parameter bx is the contribution to the effective potential from the low energy tail of 
the vector bosons 



4>/2£(3) cos 2 Owip) 2 
f 2 



= = g 3 g 10 _46 eV -4 T 5_ (n) 

7T 2 M}, 



The second term in Eq. (10) is a matter potential arising in an asymmetric background. 
The parameter b A is given by 

bA = 3 ^ (3) G F T 3 ~ 1.53 x 10- 24 eV- 2 T 3 (12) 

47T 4 



and we have defined the parameter 



*p* = jp T) . (i3) 



The term L^' T ^ contains all the asymmetries in the problem 

L< T > = 2L Pt + L Ull + L Vc - X -L n . (14) 

An asymmetry is defined as L Va = (N Ua — N-p a )/N^, where N Ua = (3/8)N 7 n Ua is the 
number density of a given neutrino specie and iV 7 = 2£(3)T 3 /vr 2 is the photon number 
density. Specifically for v T one finds 

L Vr = \{P Z -P Z \ (15) 

such that 

5P Z = P Z -P Z + P C , (16) 
where P c = (8/3) L m is a constant initial asymmetry. 



The expansion of the universe is given by the Friedmann equation H = y/8irGp/3, 
where p = ^g*T 4 . In accordance with our assumption of non-adiabatic oscillations we take 
g* = Const = 10.75 [|||. 

Since we are focusing on the evolution of small asymmetries it is an advantage to redefine 
the polarisation vectors (P and P) as 

P± = P i ±P l (i = x,y,z). (17) 

in order to separate small (P~) and large (P + ) quantities. This greatly increases the 
numerical accuracy. 

Numerical results — We have numerically studied the system of equations at temper- 
atures below 100 MeV. At high temperatures the v T will be in thermal equilibrium via its 
weak interactions. The v s can only be brought into thermal equilibrium by oscillations. 
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However, at high temperatures the damping is very strong and oscillations therefore to- 
tally washed out. Initially the abundance of v s will therefore be practically zero. The 
initial lepton asymmetry is plausible comparable in magnitude to the baryon asymmetry 
(L m ~ 10 -10 — 10 -9 (however, this need not be the case). The evolution of the leptonic 
asymmetry depends upon the sign of 5m 2 . If 5m 2 < the system encounters a MSW reso- 
nance where the effective matter mixing angle becomes maximal at a temperature (see e.g. 

H) 

T res = 16.0 (cos 29 ■ (|5m 2 |/eV 2 )) 1/6 MeV. (18) 

We find solutions very similar to those discussed by Enqvist et al.|14]. We also find that in 
large regions of parameter space the sign of the final asymmetry is extremely sensitive to 
small variations in mixing parameters and initial asymmetry. It is interesting to investigate 
if the observed sign indeterminacy is a real chaotic phenomenon in the mathematical sense 
that all information about the initial condition is lost during the process of asymmetry 
amplification or if the system is just "sensitive" . 

In order to investigate these dynamical aspects of the system, we have numerically 
calculated the Lyapunov characteristic exponents, Aj (LCE's). The LCE's measure the 
exponential divergence and convergence of nearby orbits in phase space. A dissipative 
system with a positive LCE is defined as chaotic since it ultimately losses all information 
about the initial condition. A dissipative system contracts phase space volume in time 
< 0)- The set of QRE's is a dissipative dynamical system because of the expansion 
of the universe, so that the LCE's is an ideal way to study the properties of the system. 

Xi > is not only a necessary condition for chaos, it also quantifies the degree of chaos 
in the system. Its magnitude is a measure of how long it takes for an initial uncertainty do 
develop into total chaos. The information loss / related to a postive LCE is basicly after a 
time t equal to / = At and is typically measured in bits. 

We have numerically calculated the LCE's for the QRE's ©. The calculation has been 
based on a method by Benettin [24]] . An orbit in 6-dimensional phase-space is defined by 
integrating the non-linear system of 6 QRE's. Surrounding an initial point on this reference 
orbit is fixed a sphere defined by 6 principal axes. The surface of this sphere now defines the 
different initial conditions in the problem. When integrating the non-linear system forward 
in time the sphere follows the point on the reference orbit that it is fixed to. While doing 
this it becomes deformed to a 6-ellipsoid because of the divergence or convergence of close 
orbits in phase-space. We define the LCE's from the change in the magnitude of the lengths 
\vi\ of the spheres' principal axes 

A,= lim ylog 2 j^i, (19) 



The evolution of the sphere's principal axes is defined by a linearised version of the non 
liniar set of equations. In each integration point the magnitude of the sphere's principal 
axes are compared to the magnitude at the previous integration point. Averaged over time 
this produces the time averaged LCE's [23|. 



The system of 6 QRE's contains 6 different LCE's. Until the system hits the resonance 
all 6 exponents are negative and no information lost. Once the system encounters the reso- 
nance, the leptonic asymmetry begins to oscillate rapidly from positive to negative values. 
At this point the dynamics of the system changes; it now has 1 positive and 5 negative 
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Figure 1: Upper panel: The evolution of L u with temperature for the parameters 5m? = 
— 10~ 3 eV 2 , sin 2 29 = 10~ 5 and L™ = 10~ 10 . Lower panel: The total loss of information in 
bits for the same model. 



exponents, so information is lost. Fig. 1 has been produced form the initial conditions 
(5m 2 , sin 2 29, P c )=(— 10~ 3 eV 2 , 1CP 5 , 1CP 10 ). The upper panel shows the asymmetry as a 
function of temperature inside the resonance, whereas the lower displays the loss of infor- 
mation through the resonance. This information loss has been averaged over 100 different 
initial conditions which lie very close in phase space. This eliminates noise and does not in 
any way distort the conclusions. 

The LCE's in general depend upon the initial conditions. Only if an ergodic measure 
exists for the attractor, will almost all initial conditions in the basin of attraction result in 
the same LCE's [27]. Often, the structure of the basin of attraction is so complicated that 
even a small variation in initial conditions puts one in the basin of attraction of another 
attractor, characterized by a different set of LCE's [25]. 

In Fig. 1 it seems reasonable that the situation with a negative final asymmetry can 
not be characterized by the same attractor as the situation with a positive final asymmetry. 
Most likely a different type of attractor also appears when the system crosses the resonance, 
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namely a chaotic attractor p6|]. Since different attractors in general are characterized by 
different LCE's, one must average the LCE's over different initial conditions to probe the 
spatial structure of the attractor. 

Prom Fig. 1 it is noticed that as the asymmetry initially grows and starts oscillating 
rapidly, information is lost. After the initial violent growth the lepton number oscillates 
more regularly, and information is no longer lost. Very suddenly the system begins then 
to oscillate with much smaller amplitude around one of the two attractors. This can be 
considered a chaotic transition where, as displayed in the lower part of Fig. 1, information is 
rapidly lost. The system goes through this alternating behaviour of "ordinary oscillations" 
followed by chaotic transitions all the way trough the resonance^. However, ultimately the 
system is driven out of the chaotic region by the expansion of the universe such that only 
a limited amount of information is lost. 

The system is not truly chaotic in the sense that it has a positive time averaged Lyapunov 
exponent at all times. However, it goes trough the chaotic region and will therefore be 
sensitive to variations in initial conditions. This sensitivity depends on the amount of 
information lost, which again depends critically on the value of the mixing parameters. 
Fig. 2 has been produced from the initial conditions (5m 2 , sin 2 26, P c )={— 10~ 3 eV 2 , 10 -7 , 
1CP 10 ) and illustrates nicely that the loss of information is smaller for smaller mixing angles. 

In general we find that the sensitivity grows with larger mixing angles and smaller 
mass differences^. This is clearly seen in Fig. 3 where we have calculated the total loss of 
information (in bits) for different values of 5m 2 and sin 2 26. 

Discussion — We have performed a careful study of the dynamical properties of the 
quantum rate equations for active-sterile neutrino oscillations in the early universe. Nu- 
merically we have confirmed the results of Ref. [|l4|], that there are regions of parameter 
space where the final sign of the lepton asymmetry is extremely sensitive to initial condi- 
tions. However, even though the systems goes through chaotic transitions, it is never truly 
chaotic in a mathematical sense because the system only stays a finite time in the chaotic 
regions. 

Effectively though, the system can behave exactly like a chaotic system in some regions 
of parameter space. Thermal fluctuations in the early universe can be expected to produce 
slightly different initial asymmetries, of the order AL m ~ 10~ 18 [15]. If the average initial 
asymmetry is of the order L m ~ 10~ 10 the information loss needs only be of the order 25 
bits in order to produce a completely arbitrary final sign. We have indicated this region of 
(5m 2 , sin 2 2#)-space in Fig. 3 by a light shading. 

It should be noted that there are reasons for believing that much larger initial differences 
in L could be present. They could for instance have been produced during the electroweak 
phase transition (there could be cases where AL m ~ L m ). In that case, large regions in the 
mixing parameter space leads to arbitrary sign, essentially the asymmetry only has to go 
through a few oscillations. 

The outcome of BBN depends strongly on the sign of a possible lepton asymmetry, 
if the asymmetry is in the electron neutrino sector. The reason is that electron neutrinos 



6 For larger 6m 2 the system only goes through chaotic transitions at the beginning of the resonance. 
Afterwards the asymmetry oscillates with slowly increasing asymmetry until the end of the resonance. This 
is the behaviour seen in Fig. 1 of Ref. JlJ]. However, information is only lost when the chaotic transitions 
are present. 

7 This conclusion is supported by Fig. 2 of Ref. |TJ|]. 
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Figure 2: Upper panel: The evolution of L u with temperature for the parameters 5m? = 
— 10~ 3 eV 2 , sin 2 29 = 10 -7 and = 10~ 10 . Lower panel: The total loss of information in 
bits for the same model. 

participate directly in the weak reactions that interconvert neutrons and protons. A positive 
L leads to an decrease in Yp and a negative L to an increase. So in principle different regions 
could have different light element abundances. The uncertainty in Yp is expected to be of 
the order 10~ 2 . The formation of domains with different lepton number can in principle 
also lead to visible effects in the CMB. The extra energy density produced by repopulation 
is in most cases very small (which is exactly why our non-adiabatic approximation applies) , 
so the oscillations are not detectable in this way [pfifl . However, the recombination history 
depends sensitively on Yp, so the CMB anisotropies are still altered to some extent by the 
oscillations. 

Finally, we caution that our calculations have been performed using the quantum rate 
equations (QRE's) instead of the full set of quantum kinetic equations (QKE's). The 
interaction between different momentum states has been neglected in the QRE's and it is 
likely that a full solution of the QKE's would find a somewhat smaller degree of chaos. 
However, it also seems well-established that the asymmetry does indeed oscillate also when 
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Figure 3: The loss of information, calculated for different values of 5m 2 and sin 2 29. In all 
cases we have assumed an initial asymmetry of L™ = 10 -10 . 



the QKE's are solved |15, |jq| , and so one may expect our analysis to qualitatively also apply 
to the QKE's. 
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